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ABSTRACT 

The relationship between galaxy star formation rates (SFR) and stellar masses (M*) is re-examined 
using a mass-selected sample of ~62,000 star-forming galaxies at z < 1.3 in the COSMOS 2-deg 2 field. 

Using new far-infrared photometry from Herschel- PACS and SPIRE and Spitzer-MIPS 24 /im, along 
with derived infrared luminosities from the NRK method based on galaxies’ locations in the restframe 
color-color diagram ( NUV — r) vs. (j— K), we are able to more accurately determine total SFRs for 
our complete sample. At all redshifts, the relationship between median SFR and M* follows a power- 
law at low stellar masses, and flattens to nearly constant SFR at high stellar masses. We describe a 
new parameterization that provides the best fit to the main sequence and characterizes the low mass 
power-law slope, turnover mass, and overall scaling. The turnover in the main sequence occurs at a 
characteristic mass of about M 0 ~ 1O 1O M 0 at all redshifts. The low mass power-law slope ranges 
from 0.9-1.3 and the overall scaling rises in SFR as a function of (1 + z) 4 ' 12±0,10 . A broken power-law 
fit below and above the turnover mass gives relationships of SFR oc jVf°- 88±0 - 06 below the turnover 
mass and SFR oc M° ,27±0,04 above the turnover mass. Galaxies more massive than M* > 10 10 M 0 
have on average, a much lower specific star formation rate (sSFR) than would be expected by simply 
extrapolating the traditional linear fit to the main sequence found for less massive galaxies. 

Subject headings: galaxies: evolution - galaxies: high-redshift 


1. INTRODUCTION 

Over the last decade, a tight correlation between a 
galaxy’s star formation rate (SFR) and its stellar mass 
(M*) has been discovered (Noeske et al. 2007; Daddi 
et al. 2007; Elbaz et al. 2007). Commonly referred to as 
the galaxy “main-sequence” (MS) of star formation, this 
relationship has important implications for the physical 
nature of star formation in galaxies. The MS is generally 
described as a single power law of the form SFR oc Mj 
with /3 = 0.7-1.0 and the normalization of the MS evolv¬ 
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ing to higher values at increasing redshift (Noeske et al. 
2007). 

A common interpretation of the existence and tight¬ 
ness of the main-sequence is that the majority of star¬ 
forming galaxies are powered by similar quasi-steady pro¬ 
cesses, with only a small fraction of galaxies undergoing 
more chaotic processes such as major merger events that 
might be expected to produce strong bursts of star for¬ 
mation (e.g. Elbaz et al. 2011; Rodighiero et al. 2011; 
Sargent et al. 2012). These starburst galaxies are gener¬ 
ally thought to lie significantly above the MS and repre¬ 
sent a minority of galaxies. 

A key uncertainty in measuring galaxy SFRs is the 
effect of dust obscuration. The most direct method of 
determining dust obscuration is from observations in the 
far-infrared, where the absorbed starlight is thermally 
reradiated. In the absence of far-infrared data, various 
extrapolations from shorter wavelength have been used 
to study the main sequence, such as using emission lines 
combined with reddening corrections to infer the dust- 
corrected SFR (Brusa et al. 2010; Sobral et al. 2012; Za¬ 
hid et al. 2012; Kasliino et al. 2013) or measuring the 
UV or optical emission from young massive stars and 
correcting for the radiation lost to dust obscuration (Lee 
et al. 2011; Rodighiero et al. 2011; Steinhardt et al. 2014). 
Observations in the mid-infrared (e.g. 24 pm) have been 
used to estimate far-infrared luminosities (e.g. Noeske 
et al. 2007; Daddi et al. 2007; Elbaz et al. 2007), al¬ 
though the accuracy of these estimates decreases at high 
redshifts and bright infrared luminosities (e.g. Papovich 
et al. 2007; Lee et al. 2010; Elbaz et al. 2011). Stud¬ 
ies of radio emission take advantage of the well-known 
radio-FIR correlation (Helou et al. 1985; Condon 1992; 
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Yun et al. 2001) to estimate the infrared luminosities, but 
many of these studies rely on stacking to overcome high 
sensitivity limits (Dunne et al. 2009; Pannella et al. 2009; 
Karim et al. 2011). The consensus from these studies is 
that the MS follows a single power law SFR oc M* , with 
the slope generally between j3 = 0.7-1.0 and the normal¬ 
ization varying based on the study’s redshift, SFR in¬ 
dicator, sample selection, and IMF (for a summary, see 
Speagle et al. 2014). 

However, a few studies have found indications of a 
more complex main-sequence relationship. Some stud¬ 
ies suggest that the MS slope varies with stellar mass 
so that a single power-law cannot explain the MS and 
a stellar mass-dependent slope is a better fit (Karim 
et al. 2011; Whitaker et al. 2012; Magnelli et al. 2014). 
Recent studies based on far-infrared selected samples 
from Herschel show that far-infrared selected galaxies lie 
mostly above the MS, with a much shallower slope in 
\og(SFR)/log(M*) (Lee et al. 2013; Oteo et al. 2013a,b; 
Lemaux et al. 2013). However, this discrepancy is due 
to the flux limited selection of far-infrared samples that 
introduce a SFR-based selection bias, as compared to 
studies based on stellar mass-selected galaxy samples. 
This has been demonstrated by stacking analyses that 
explore the far-infrared emission as a function of stellar 
mass and find generally good agreement between dust- 
corrected UV-derived SFRs and Herschel -derived SFRs 
(e.g. Rodighiero et al. 2014). 

Stacking is a commonly used technique to measure low- 
level emission from galaxies that would be undetected 
individually. Stacking analyses require a number of as¬ 
sumptions and can miss vital information about indi¬ 
vidual galaxies and their distributions. Unless the par¬ 
ent population is identical (a key assumption in stack¬ 
ing), interpretation of stacking results can be difficult be¬ 
cause the underlying distribution is unknown (although 
see Schreiber et al. 2014, for a possible method to de¬ 
termine the underlying distribution). In addition, these 
stacking analyses do not explain why the dust-corrected 
SFRs cannot accurately recover the SFRs seen in high 
luminosity galaxies, which have an elevated contribution 
to the integrated build up of stellar mass in the universe. 

Direct Herschel FIR measurements remain a unique 
tool to properly estimate the ongoing star formation rate 
in the most active dusty galaxies. Analysis of the rest- 
frame UV emission in dusty galaxies suggests that apply¬ 
ing the nominal attenuation laws (e.g. Meurer et al. 1999; 
Calzetti et al. 2000) will dramatically underestimate to¬ 
tal star formation rate in galaxies exceeding ~ 5OM 0 /yr 
(Smail et al. 2004; Casey et al. 2014; Rodighiero et al. 
2014). 

In the following paper, we attempt to address these is¬ 
sues by using a dust-corrected SFR indicator that is accu¬ 
rate for galaxies at all luminosities. By analyzing a large 
sample of individual galaxies, we do not lose information 
about the distribution of sources from stacking and can 
re-examine the shape of the star-forming MS in a stellar 
mass-selected sample. The data are described in Sec¬ 
tion 2 and SFRs computed by several different methods 
are measured and compared in Section 3. In Section 4 
we analyze our mass-selected sample of galaxies in the 
SFR/M * plane and find the best fits to the data. The 
implications of the main-sequence are discussed in Sec¬ 


tion 5 and we list our conclusions in Section 6. When cal¬ 
culating rest-frame quantities, we use a cosmology with 
f l m = 0.28, A = 0.72, and H 0 = 70 km s -1 Mpc -1 
(Hinshaw et al. 2013). A Chabrier (2003) Initial Mass 
Function (IMF) truncated at 0.1 and 100 M 0 is used 
when deriving SFRs and stellar masses. 

2. DATA 

Our analysis of the MS is made possible by the large 
area multi-wavelength coverage of the COSMOS field, a 
2 deg 2 area of the sky with observations from the ul¬ 
traviolet through far-infrared and radio (Scoville et al. 
2007). We construct a mass-complete sample of galax¬ 
ies with K s < 24 from the deep AT s -band catalog of II- 
bert et al. (2013), based on data from the first Ultra- 
VISTA DR1 data release, covering -^75% of the COS¬ 
MOS field (McCracken et al. 2012). 20 bands of opti¬ 
cal and near-infrared photometry were extracted using 
matched apertures in dual-image mode from the vari¬ 
ous available COSMOS images and was combined with 
GALEX magnitudes from the multi-wavelength catalog 
of Capak et al. (2007). We use updated Spitzer IRAC 
photometry from the Spitzer Large Area Survey with 
Hyper-Suprime-CAM (SPLASH, Capak et al. in prep). 
We cross-match this catalog with the Spitzer MIPS 24 
p,m catalog of Le Floc’h et al. (2009) and the Herschel 
catalog of Lee et al. (2013) using a matching radius of 
2 ". 

2.1. Source Selection 

We interpolate the 90% stellar mass completeness lim¬ 
its from Ilbert et al. (2013) to determine approximate 
mass completeness thresholds at all redshifts, and se¬ 
lect only galaxies with stellar masses above their red- 
shift dependent mass completeness limit. When study¬ 
ing star-forming galaxy populations, we separate “star- 
forming” and “quiescent” galaxies using a two-color se¬ 
lection technique: NUV— r + versus r + — J (as described 
in Ilbert et al. 2013). Specifically, galaxies with absolute 
magnitude colors Mnuv — M r > 3 (M r — Mj) + 1 and 
Mnuv ~ M r > 3.1 are considered “quiescent” (~ 15% of 
the sample) while the remaining galaxies are considered 
actively star-forming galaxies. 

Star-forming galaxies that also contain luminous AGN 
are a concern because the luminosity from the AGN is 
extremely difficult to separate from emission from star- 
formation, and thus these sources may have erroneously 
high SFRs (although this concern is lessened for FIR 
sources because AGN generally heat dust to tempera¬ 
tures too hot to radiate in the far-infrared). On the 
other hand, many of the galaxies that host AGN also 
contain significant star-formation, and removing these 
sources introduces a bias to our study. We find that the 
overall results of our study are not significantly affected 
by either the inclusion or exclusion of these sources, so 
we do not remove galaxies that have been detected in 
the X-ray (~ 0.5% of sample) by XMM-Newton (Brusa 
et al. 2010) or Chandra (Civano et al. 2011), or that have 
IRAC power-law colors (Donley et al. 2008) that suggest 
AGN activity (~ 25% of sample). 

2.2. Infrared Data 

The Herschel-selected sample of galaxies is described 
in detail in Lee et al. (2013, hereafter L13) and is briefly 
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summarized here. L13 use Spitzer 24 /tin and VLA 
1.4 GHz priors to find 4,218 sources in COSMOS that 
were each detected in at least two of the five available 
Herschel PACS (100 /tm or 160 /tm) and SPIRE (250 
/tm, 350 /inr, or 500 /tm) bands. These sources span 
log (Air/A q ) = 9.4-13.6 and 2 = 0.02-3.54. Dust prop¬ 
erties of each source (e.g. Air, Adust, A/dust) were mea¬ 
sured by fitting the full infrared photometry to a coupled 
modified blackbody plus mid-infrared power law using 
the prescription given in Casey (2012) and assuming an 
opacity model where r = 1 at 200 /tm. 

There is a population of galaxies that are classified as 
“quiescent” from their NUV—r + versus r + — </ colors, 
but have been detected in the infrared by Herschel or 
Spitzer , suggesting that these galaxies are actually un¬ 
dergoing a significant amount of star-formation (~ 7% 
of the “quiescent” population). These galaxies are more 
consistent with being very dusty objects that have ex¬ 
tremely red colors due to obscuration, not lack of active 
star-formation, so we include these galaxies in our sample 
of star-forming galaxies (see Section 4.3). 

2.3. Photometric Redshifts and Physical Parameters 

Ilbert et al. (2013) measure accurate 30-band pho¬ 
tometric redshifts of the full IQ-band COSMOS cata¬ 
log. We find a median Az/(1 + z) = 0.02 in our sam¬ 
ple of star-forming galaxies, with a catastrophic failure 
(|Az|/(l + z) > 0.15) in 5.6% of sources. In addition, 
physical parameters such as stellar mass and star forma¬ 
tion rate have been calculated by fitting the Spectral En¬ 
ergy Distributions (SEDs) to synthetic spectra generated 
using the Stellar Population Synthesis (SPS) models of 
Bruzual & Chariot (2003). We also recalculate the phys¬ 
ical parameters using different templates and extraction 
parameters, and find that our final results are not af¬ 
fected by the specific choice of template. Thus, we use 
the same set of parameters used to create the catalog in 
Ilbert et al. (2013), but with updated near-IR photome¬ 
try from SPLASH. 

3. ANALYSIS 

3.1. Star Formation Rate Calculations 

There are many methods for estimating a galaxy’s SFR 
based on observations at various wavelengths (for a re¬ 
view, see Kennicutt 1998; Murphy et al. 2011). Here we 
compare a few commonly used SFR indicators using a 
subset of COSMOS galaxies to determine how much the 
different SFR methods disagree. In all cases, we measure 
the total SFR as SFR Tot = SFR m + SFR Vv . 

3.1.1. Infrared derived SFR 

As discussed in Section 2.2, the infrared properties of 
the Herschel- selected galaxies have been measured by fit¬ 
ting the infrared SEDs to a coupled modified blackbody 
plus mid-infrared power law model (Casey 2012). It has 
been shown that measuring the Air from fitting the far- 
infrared data to libraries of SED (e.g. Chary & Elbaz 
2001; Dale & Helou 2002) gives roughly the same results 
as the modified blackbody plus power-law model (Casey 
2012; U et al. 2012; Lee et al. 2013). The infrared ob¬ 
servations give us an estimate of the obscured SFR, and 
we combine this with UV observations of the unobscured 


SFR to derive the total SFR as in Arnouts et al. (2013): 

SFR t otal = (8.6 x 10 -11 ) x (Air+2.3x rA„(2300A)) (1) 

where Air = A(8-1000/mi) and all luminosities are mea¬ 
sured in units of A 0 . For sources with SFR. > 5OM 0 /yr, 
the infrared contribution dominates the total SFR, con¬ 
tributing as much as ~ 90% of the total SFR. 

While we have excellent Herschel coverage of the full 2- 
deg 2 COSMOS field that yields 4,218 sources, the detec¬ 
tion limits of Herschel introduce a selection bias against 
all but the most luminous infrared sources. A common 
method of determining the Air of less luminous galax¬ 
ies is to use deep Spitzer 24 /tm data to estimate the 
far-infrared luminosity (e.g. Kennicutt et al. 2009; Rieke 
et al. 2009; Rujopakarn et al. 2013). COSMOS has ex¬ 
tremely deep coverage at 24 /tm and Le Floc’h et al. 
(2009) provide SFR estimates for 36,635 galaxies, which 
we use to extend our sample of infrared detected galaxies 
to more moderate luminosities. 

3.1.2. Optical & UV based SFR Indicators 

For galaxies without direct measurements from far- 
or mid-infrared wavelengths of the obscured SFR, the 
amount of radiation obscured by dust must be estimated 
indirectly. A common method for estimating total SFR 
is to fit libraries of model SEDs (that include prescrip¬ 
tions for dust obscuration) to optical & UV photometry. 
Ilbert et al. (2013) use the full optical COSMOS photom¬ 
etry and fit to a library of synthetic spectra from Bruzual 
& Chariot (2003), and estimate the total SFR for each 
of the galaxies in our sample from the best fit SEDs. 

Another method of estimating dust-corrected SFRs is 
by using rest-frame UV observations to measure the un¬ 
obscured SFR and inferring the appropriate dust correc¬ 
tion factor from observed colors. Two examples of this 
are the BzK method from Daddi et al. (2004) and the 
NRK method from Arnouts et al. (2013). BzK SFRs are 
determined by using the observed-frame R-band photom¬ 
etry to measure the rest-frame UV luminosity, and then 
estimating the extinction as E{B — V) = 0.25(1? — z + 
0.1 )ab (Daddi et al. 2007). BzK SFRs are only valid for 
redshifts 1.4 < z < 2.5, as these are the only redshifts 
where the desired portions of the SED are redshifted 
to the correct wavelengths, and we limit our selection 
to the good- sBzK with errors <51og[SFR(UV)]< 0.3 dex 
(Rodighiero et al. 2014). 

NRK SFRs are calculated by using their location in the 
rest-frame color-color diagram ( NUV — r) vs. (r — K) 
to estimate extinction. Arnouts et al. (2013) find that at 
z < 1.3, the infrared excess IRX = Air/ Lnuv in star¬ 
forming galaxies can be parameterized as a function of 
redshift and the vector NRK= 0.31 x (NUV —r) + 0.95 x 
(r — K). This allows us to estimate the Air and calcu¬ 
late total SFR using Equation 1. When measuring NRK 
SFRs, we use the small “sSFR correction” as described 
in Arnouts et al. (2013). 

3.2. Comparison of SFR indicators 

We compare commonly used SFR indicators using a 
common subset of COSMOS galaxies to determine how 
much agreement there is between the different measures 
of SFR. We have a large set of Herschel detected galax¬ 
ies from Lee et al. (2013) where, for the first time, we 
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have direct measurements of both the obscured and un¬ 
obscured SFR (from UV observations) at a wide range 
of redshifts. We compare the other SFR indicators dis¬ 
cussed previously (24 /im, SED fits, BzK, and NRK) to 
this sample of 4,218 Herschel detected galaxies. 

Figure 1 displays the comparison of SFR T otal from the 
four different indicators discussed above to SFR.T 0ta \ as 
measured by Herschel. Density contours show the loca¬ 
tion and concentration of the majority of the sources, 
with outliers shown in gray circles. Median values 
in 20 equally populated bins of SF Rectal ,Herschel are 
over-plotted to show average trends. To determine the 
strength of the correlation between each SFR indicator 
and S'FRxotai, Herschel, we measure the Pearson correla¬ 
tion coefficient (p) and provide these values at the top 
of each sub-panel. The Pearson correlation coefficient 
can vary between +1 and -1, with +1 indicating total 
positive correlation, 0 indicating no correlation, and -1 
indicating total negative correlation. We also measure 
the median difference between each SFR indicator and 
Herschel SFR (< A log(SFR) >) and list these values at 
the top of each sub-panel. 

The 24 /im-determined SFR correlates with the Her¬ 
schel SFR very well (p 24 = 0.88, < Alog(SFR 2 4 ) > 
= 0.12), except at the highest IR luminosities. This trend 
has been previously explored in many studies which find 
that at moderate redshifts and IR luminosities, 24 /im 
observations are a good proxy for im , but at high red¬ 
shifts and infrared luminosities, the 24 //in estimates tend 
to overpredict the true Air, possibly due to redshifting of 
the observed 24 //m-band to wavelengths contaminated 
by PAH features (e.g. Papovich et al. 2007; Lee et al. 
2010; Elbaz et al. 2011). As we are using the 24 /zm 
SFRs to fill in the low and moderate luminosity galaxies 
that are not detected with Herschel , this discrepancy is 
not a major issue for our work. 

The NRK SFRs also show strong correlation with 
the Herschel- derived SFRs ( Pnrk = 0.79, < 
A log(SFRNRK) > = 0.17). This is not completely un¬ 
expected since the NRK method was developed using 24 
//m-derived SFRs as a baseline, but the NRK measured 
SFRs match very well with those derived from Herschel. 
Like with SFR 24 , the correlation shows signs of break¬ 
ing down at the highest SFRs, but as long as the NRK is 
used mainly for low SFR galaxies, it provides a reliable 
estimate of the SFR. 

By contrast, the agreement between SFRsed and 
SFRt otal is quite poor, showing much weaker correla¬ 
tion between the two indicators (psed = 0.56). The 
tightness of the correlation is also much broader (< 
Alog(SFRsED ) > = 0.43), even at low SFRs where 
the median points lie closer to the unity line. Again, 
at high SFRs the median points show a clear deviation 
from unity. Wuyts et al. (2011) are able to find a bet¬ 
ter match between SFRsed and SFR 24 if they tune key 
parameters of the SED fit, such as r m j„, the e-folding 
time of the exponentially declining star formation his¬ 
tory. The exact tuning needed varies based on several 
other assumptions in the SED fitting procedure, such 
as different stellar population synthesis codes, and even 
when the tuning is done, the computed SFR^ee still sys¬ 
tematically underestimates the true SFR for a significant 
fraction of sources. 

Finally, the comparison between SFRbzK and 


S ERxotai shows essentially no correlation (PBzK = 
—0.11). It should be noted that the redshift range of 
the BzK indicator (1.4 < z < 2.5) limits us to a small 
sample size containing only the brightest galaxies. As 
seen in Figure 2, this selection limits our comparison 
to galaxies at SFRs where all indicators begin to devi¬ 
ate significantly from S'Fi?Totai- Stacking analyses sug¬ 
gest a stronger correlation between average SFRbzK and 
average SFRnerschei at fainter luminosities (Rodighiero 
et al. 2014), but the tightness of the distribution is not 
well determined. The BzK galaxies that are Herschel 
detected show no correlation between the SFR derived 
from the BzK method and from Herschel measurements. 

3.2.1. Selection Effects of SFR Indicators 

The comparisons of the various SFR indicators shown 
in Figure 1 span different dynamic ranges in SFR, mostly 
due to the redshift limitations of the NRK and BzK 
indicators. In Figure 2, we plot the typical difference 
between SFR indicators and SFRt otal as a function of 
S ERxotai • We see that all of the SFR indicators pro¬ 
vide poor estimates of the infrared measured SFRt otal 
above log(SFR) > 1.5 (~ 3OM 0 /yr). These common 
SFR indicators fail to accurately estimate the true SFR 
of luminous infrared galaxies. This highlights the need 
for direct infrared observations to accurately measure the 
SFR of highly star forming galaxies. 

Throughout this analysis we have assumed the 
Herschel-detemained SFR is the most accurate because it 
directly probes far-infrared wavelengths, where the bulk 
of the re-radiated radiation from dust is emitted. How¬ 
ever, it is possible that the high detection threshold of 
Herschel limits us to a biased sample that does not accu¬ 
rately reflect the emission properties of lower luminosity 
galaxies. To test this possibility, we re-run our analyses 
using the much deeper sample of Spitzer 24 //m-detected 
galaxies (which showed excellent agreement with Her¬ 
schel SFRs) as the comparison sample. We find very 
similar results as the Herschel comparison, with NRK 
providing both the strongest correlation and the tightest 
distribution. 

3.2.2. A Ladder of SFR Indicators 

While all three non-infrared based SFR indicators fail 
to accurately estimate the SFR in high luminosity galax¬ 
ies, the NRK method provides the most accurate and 
consistent estimates across the full dynamical range of 
Herschel SFRs. At high SFRs ( SFR > 3OM 0 /yr), 70% 
of our sample is directly detected in the infrared by either 
Herschel or Spitzer. Thus, we can study the full popula¬ 
tion of star forming galaxies by constructing a “ladder” 
of SFR indicators (as in Wuyts et al. 2011) based on 
the Herschel , Spitzer 24 p m, and NRK SFR indicators. 
All sources have SFR T otal calculated using Equation 1, 
with different methods of determining Air. For sources 
detected by Herschel , we measure Air from fitting the 
far-infrared photometry to the Casey (2012) greybody 
plus power-law models. We use the Air estimated from 
24 //nr (Le Floc’h et al. 2009) for sources that are not 
detected by Herschel but are detected at Spitzer 24 //nr. 
And for the remaining sources, we estimate the Air us¬ 
ing the NRK-derived IRX (as discussed in Section 3.1.2). 
Although we include NRK-derived IRX for all galax¬ 
ies above our mass-conrpleteness limits, the method has 
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only been well-calibrated for M* > 1 O 93 M 0 . Infrared 
stacking suggests that any systematic offsets should be 
small, but when calculating main sequence relationships 
we only include galaxies with M* > 10 9 3 Mg. 

The relative fraction of sources with SFRs measured 
from each indicator is plotted in Figure 3 as a function 
of both stellar mass and redshift. Below M* < 10 9 , 5 Mq, 
SFRs are almost all determined from NRK, but at higher 
stellar masses, the fraction of sources with direct infrared 
measurements increases until about 25% (60%) of sources 
at M > 10 10 5 Mq have SFRs determined from Herschel 
(Spitzer 24 ;tm). Because of the redshift limitations of 
the NRK method (see Arnouts et al. 2013) and the larger 
errors associated with the SFRsw and SFR^ z k indi¬ 
cators, we restrict the rest of our analysis to redsliifts 
0 < z < 1.3, where we can more accurately measure the 
SFR of our full sample. 

4. SHAPE OF THE MAIN SEQUENCE OF STAR 
FORMATION 

With reliable and consistent SFR estimates for a large, 
mass-complete sample of galaxies in COSMOS, we exam¬ 
ine the star-forming main sequence for a large, unbiased 
sample of 62,521 galaxies. Figure 4 displays the stellar 
mass and SFR of our full sample, split into four redshift 
bins spanning 0.2 < z < 1.3. Black contours display 
the density of sources at each location in SFR and M* 
parameter space, and colored bars represent the median 
SFR in stellar mass bins of width Alog(M*) = 0.3, with 
vertical error bars displaying the standard deviation of 
the SFRs in that bin. These same bins are used to create 
the fractional histograms plotted on the side of each red¬ 
shift bin, which display the distribution of SFRs within 
each mass bin with the corresponding color. The derived 
main sequence relationships from star-forming galaxies 
in Karim et al. (2011) and Whitaker et al. (2014) are 
also plotted for comparison. 

Figure 4 shows that the galaxies in our sample do 
not follow a simple linear main sequence relationship 
between log(SFR) and log(AR) (or a single power-law 
relationship between SFR and M*). Instead, the me¬ 
dian SFR relationship appears to flatten at masses above 
M* ~ 10 1 o Mq. This can be seen in the histograms, 
which show that the peak of each SFR distribution in¬ 
creases with increasing M* at low masses, but at high 
masses, the histogram peaks all lie at approximately the 
same SFR. The standard deviation of the SFR in each 
stellar mass bin remains mostly constant at all masses 
and at all redshifts, with a ~ 0.36 dex in all bins. The 
shape of this relationship appears roughly constant with 
redshift, with the entire relationship increasing to higher 
SFRs at higher redshifts. 

4.1. Parameterizing the Star-Forming Main Sequence 

From Figure 4, it is clear that a single power-law does 
not accurately describe the relationship between stellar 
mass and star formation rate. We split our sample of 
star-forming galaxies into 6 equally populated redshift 
bins, each of which are then split into 30 equally popu¬ 
lated stellar mass bins (with ~ 350 sources in each bin), 
and calculate the median SFR in each bin. We limit our 
sample to stellar masses above a conservative mass limit 
(see Table 1) to ensure that we are not affected by sys- 
tematics. The specific number of redshift bins does not 


affect the following results, although we must balance be¬ 
tween having redshift bins that are too wide and combine 
galaxy samples at different epochs with having redshift 
bins that are too narrow and are affected by small num¬ 
ber statistics. The same is true for the number of stellar 
mass bins, although having at least 30 bins is prefer¬ 
able for accurately determining the goodness of fit to the 
models. The median SFRs for every bin are plotted in 
Figure 5, colored by redshift, with bootstrapped errors 
on the median represented by vertical bars. 

Using the MPFIT package implemented in IDL (Mark- 
wardt 2009), we fit the median log(M^) and log(SFR) in 
each redshift bin with many models including linear, 2 nd 
order polynomial, and broken linear, and find the best 
fit is provided by the following model: 


S = Sq — log 


1 + 



( 2 ) 


where S = log(SFR) and M = log^M^/Mo). We choose 
this model because ( 1 ) at all redshifts, it provides the 
best reduced y 2 fit to the data, and ( 2 ) unlike polyno¬ 
mial fits, the parameters of the model allow us to quan¬ 
tify the interesting characteristics of the relation between 
stellar mass and SFR: 7 , the power-law slope at low stel¬ 
lar masses, Mo, the turnover mass (in log(M* / M@)) , and 
So, the maximum value of S (or the maximum value of 
log(SFR)) that the function asymptotically approaches 
at high stellar mass. The best-fit parameters for each 
redshift bin are listed in Table 1. 


4.2. Evolution of Model Parameters 

In the top panels of Figure 6 , we plot the evolution 
of So, Mo, and 7 as functions of log( 1 + z). The bot¬ 
tom panels of Figure 6 examine the covariance between 
these parameters by displaying the 95% confidence error 
ellipses. 

We see clear and strong evolution of S 0 with redshift, 
and the best fit line suggests an evolution of S 0 oc (4.12 ± 
0 . 10 )xZog(l+ 2 ), or equivalently, SFR Q oc (l+z) 412±01 °. 
The covariance between S 0 and Mo (Figure 6 D) and 
between So and 7 (Figure 6 E) is relatively minor, so we 
infer that the evolution in So is true evolution and not 
due to variation in the other parameters. 

Both Mo and 7 show some evidence for weak evolu¬ 
tion to more massive Mo and steeper 7 with redshift, 
although much of the perceived evolution may be due 
to the covariance seen in Figure 6 F. The best (linear) 
fit to the evolution in the turnover mass is given by 
Mo oc (1.41±0.20) x log(l + z). We test the possible red¬ 
shift evolution of turnover mass by calculating where the 
data deviates by 0.2 dex from a single power-law fit to 
the low mass data and find similar evolution, suggesting 
that the turnover mass does indeed change with cosmic 
time. The low-mass power-law slope, 7 , has a best-fit line 
that suggests evolution of 7 oc (1.17 ± 0.13) x log( 1 + z). 
Redshift evolution in 7 to steeper slopes at earlier cosmic 
times would suggest that the SFR in the lowest stellar 
mass galaxies does not increase as much as in more mas¬ 
sive systems. 

4.3. Separating Quiescent Galaxies 
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We have described the main sequence relationship be¬ 
tween SFR and M* for star-forming galaxies. How¬ 
ever, possible misclassification of galaxies as either “star¬ 
forming” or “quiescent” could drastically affect the 
trends we observe. 

As described in Section 4, we remove galaxies that are 
considered quiescent from our sample using the selection 
Mnuv — M r > 3 (M r — Mj ) ± 1 and Mnuv ~ M r > 3.1 
(Ilbert et al. 2013). This selection is shown in Figure 7, 
with the full mass-selected sample of COSMOS galaxies 
at 0.2 < z < 1.3 generally separated into two distinct 
“star-forming” and “quiescent” regions. Improper clas¬ 
sification of the “in-between” galaxies that are not obvi¬ 
ously star-forming or quiescent could lead to changes in 
the main sequence shape. To test this possibility, we shift 
the entire separating line (both horizontal and diagonal 
segments) between quiescent and star-forming galaxies 
by ±0.4 mag in Mnuv — M r , and in either case there is 
no appreciable change to the main-sequence. 

Galaxies detected in the infrared by Herschel or Spitzer 
24 pm are highlighted in Figure 7, and while the majority 
fall on the star-forming sequence, we see a number of ob¬ 
jects that lie in the quiescent region. This population of 
infrared-detected quiescent (IR-Q) galaxies is relatively 
small, with only ~ 7% of the galaxies classified as quies¬ 
cent having detectable infrared emission, but these mis- 
classified galaxies are predominantly found at high stellar 
mass. The fraction of quiescent galaxies detected in the 
infrared increases rapidly from < 1% at M* ~ 1O 9 ' 5 M 0 
to 15-20% at M* > 10 n Mg, and this trend holds at all 
redshifts. These massive galaxies could heavily influence 
the shape of the main-sequence we observe, so it is vital 
to understand what is driving their infrared emission. 

There could be several reasons why galaxies with qui¬ 
escent colors have significant emission in the infrared, in¬ 
cluding (i) improper classification of star-forming galax¬ 
ies possibly due to extreme dust obscuration, (ii) el¬ 
evated infrared luminosity from an AGN, (iii) inaccu¬ 
rate absolute magnitudes due to catastrophic failures 
in photo-z’s or low signal-to-noise photometry, or (iv) 
“post-starburst” infrared glow due to dust heating from 
young stars (that is not related to the instantaneous star 
formation). 

AGN typically heat dust to very hot temperatures, so 
we expect any AGN contribution to infrared radiation 
to be predominantly in near- and mid-IR wavelengths, 
while far-infrared emission is likely due to star formation 
alone. Only about 10% of the IR-Q galaxies have been 
detected by Herschel , and the rest are 24 pm-only detec¬ 
tions, where AGN may heavily influence the emission. 
However, only ~ 1 -5% of the IR-Q galaxies are detected 
in the X-ray by Chandra , and only ~ 2-8% of the IR-Q 
galaxies have IRAC power-law colors indicative of AGN, 
with significant overlap in those two populations, and 
the percentages are even lower when looking only at the 
24 pm-only sources. This suggests that radiation from 
an AGN is not fueling the infrared emission. The aver¬ 
age SFR of the IR-Q galaxies with AGN is 0.2-0.5 dex 
higher than the average SFR of all the IR-Q galaxies, 
so the presence of AGN in these galaxies is likely just 
a reflection of the well-studied trend that AGN fraction 
increases with SFR or Lir (e.g. Kartaltepe et al. 2010). 

The rather high SFRs of the IR-Q galaxies suggest that 
they are indeed driven by star-formation, and have been 


misclassified as quiescent. Man et al. (in prep) stack the 
infrared emission from quiescent galaxies and find up¬ 
per limits of SFRm < 0.1-1 Mg/yr. The SFRs of the 
IR-Q galaxies in our sample tend to lie below the main- 
sequence, but are all at least x2-3 higher than the upper 
limit from Man et al., which suggests that they are indeed 
still actively star-forming. In addition, the infrared emis¬ 
sion from IR-Qs is brighter than expected from a “post- 
starburst” glow (Hayward et al. 2014). It is unlikely 
that catastrophic photo-z errors or low signal-to-noise 
photometry are causing these misidentifications, as the 
sources have excellent photometry and well-constrained 
photometric redshifts (only 0.1% of the IR-Q galaxies 
have <taz/(i+z) > 0.15). Thus, the likely explanation 
for these sources is that they are actively star-forming 
galaxies that have been misclassified as “quiescent”, and 
we include them in our analysis of the main-sequence. 
We note, however, that the shape of the main-sequence 
does not change significantly based on the inclusion or 
exclusion of these sources. 

5. DISCUSSION 

We see that the slope of the main sequence relationship 
between SFR and M* changes with stellar mass. While 
most previous studies found a constant main sequence 
slope (for a summary see Speagle et al. 2014), some re¬ 
cent studies found a curved relationship might provide a 
better fit to the data (Karim et al. 2011; Whitaker et al. 
2012). However, the mass completeness limit in both 
studies coincided with the turnover mass, leaving doubt 
as to whether the turnover was a real trend or an artifact 
of completeness. With our new COSMOS observations, 
we are able to study star-forming galaxies considerably 
less massive than the turnover mass, and we find that a 
single power-law does not provide the best description of 
the star-forming main-sequence. 

5.1. The Turnover in the Star-Forming Main Sequence 

The relationship between SFR and M* varies with stel¬ 
lar mass, with two distinct regions below and above 
the characteristic turnover mass, A4o- What causes 
this change, and why does the turnover occur at about 
M* ss 10 ln Mg at all redshifts? 

5.1.1. The slope of the main-sequence 

The parameterization of the main sequence we employ 
in Section 4.1 includes a parameter 7 that we describe as 
the low stellar mass power-law slope. However, we note 
that this slope is not derived from an actual power-law fit 
to the data, but instead represents the power-law slope 
that the relationship approaches in the very low-mass 
regime, based on Equation 2. This slope is significantly 
steeper than power-law slopes commonly quoted in the 
literature, and should not be compared to slopes from 
power-law fits to data. 

For an easier comparison to the existing literature, we 
derive best-fit power-law relationships, fitting the low 
mass regime and high mass regime separately. Galax¬ 
ies less massive than the turnover mass follow a fairly 
tight power-law relationship of SFR oc M^, with /3 = 
0.88 ± 0.06. This slope is shallower than 7 because it in¬ 
cludes galaxies in the “turnover region”, where the slope 
is already starting to flatten. Galaxies more massive than 
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the turnover mass follow a drastically different relation¬ 
ship, with p = 0.27± 0.04 for M* > io 10 M Q . A galaxy’s 
specific star formation rate (SSFR = SFR/M *) can be 
interpreted as a measure of the efficiency of current star 
formation as compared to its past average star forma¬ 
tion history. The SSFR of massive galaxies is systemat¬ 
ically lower than would be expected from an extrapola¬ 
tion of low mass galaxies, suggesting that there may be 
decreased star formation efficiency in high stellar mass 
galaxies. 

5.1.2. Quenching in High Mass Galaxies 

The turnover in the main sequence to lower star for¬ 
mation efficiencies in massive galaxies suggests there is 
a fundamental change that occurs as galaxies become 
more massive, as has been predicted in some studies. 
Galaxy luminosity and mass functions, which measure 
the brightness and mass distribution of galaxies at vari¬ 
ous lookback times, show a steep, exponential decline at 
high stellar masses and high luminosities while retaining 
a remarkably consistent shape at all redshifts (e.g. Bell 
et al. 2003; Pozzetti et al. 2010; Ilbert et al. 2013). The 
lack of large, bright galaxies throughout cosmic time ar¬ 
gues for the presence of a characteristic mass above which 
a galaxy is likely to have its star formation strongly sup¬ 
pressed or quenched. 

In contrast, the dark matter halo mass function from 
semi-analytic models does not show the same exponen¬ 
tial decline, and instead has a much shallower power- 
law cutoff at much higher masses (Somerville & Primack 
1999; Benson et al. 2003), leading to a “pivot mass” 
above which the ratio of dark matter to light matter in¬ 
creases rapidly (e.g. Leauthaud et al. 2012). At low stel¬ 
lar masses, the stellar to halo mass relation (M^ oc Mf 46 ; 
Leauthaud et al. 2012) and the dark matter halo growth 
rates from IV-body simulations (Mhaio oc {Mhaio) 1 ' 1 ', 
Wechsler et al. 2002; McBride et al. 2009; Fakhouri 
& Ma 2010) suggest a main sequence relationship of 
SFR oc M 4 04 , similar to the slope seen in the main- 
sequence. The “pivot mass”, above which the stellar- 
to-halo mass relation deviates from the low stellar mass 
relationship, appears to evolve to higher M* at higher 
redshifts (Leauthaud et al. 2012), at a rate similar to the 
possible evolution seen in the main sequence turnover 
mass, Mo- In galaxies more massive than the “pivot 
mass,” the halo mass rises sharply in comparison with 
stellar mass, suggesting that while massive dark matter 
haloes appear to continue growing, the galaxies residing 
in them quench their star formation. 

Possible mechanisms for this quenching include struc¬ 
tural disruptions or galaxy mergers (Sanders & Mirabel 
1996; Hopkins et al. 2006), feedback from accretion onto 
a supermassive black hole (Springel et al. 2005), gravi¬ 
tational heating of the surrounding intracluster medium 
(Ivhochfar & Ostriker 2008), changes in the mode of gas 
accretion onto galaxies (Keres et al. 2005; Birnboim et al. 
2007; Nelson et al. 2013), or gas removal or strangulation 
in dense environments (Peng et al. 2010, 2012). 

Morphological studies may be key for understanding 
the star-formation in massive galaxies. Abramson et al. 
(2014) find that galaxy SFRs are more strongly corre¬ 
lated to disk stellar mass (as opposed to total stellar 
mass), and that SSFRdisk is approximately constant 


with mass. If this is the case, the turnover in the main- 
sequence could be simply due to growing bulges in the 
highest mass systems. However, one might expect the 
turnover to disappear (or become less severe) at high red¬ 
shifts as galaxies become more disk-dominated, but this 
is not seen in the data. Schawinski et al. (2014) find that 
disk galaxies and elliptical galaxies likely quench their 
star formation rates through different processes with very 
different timescales. A galaxy’s physical size may also 
play a role in quenching, as the surface mass density 
has been shown to correlate strongly with SSFR (Kauff- 
mann et al. 2003; Franx et al. 2008), and compact star 
forming galaxies may be on the evolutionary path toward 
quiescent galaxies (Barro et al. 2014). 

Our data suggest that galaxies with high stellar mass 
(M* > 10 10 Mq) are forming stars at a lower rate than 
would be expected from extrapolating the trends of low 
stellar mass galaxies. Finding the possible causes of this 
“quenching” of star formation is one of the key hurdles 
for understanding galaxy evolution. The existence of a 
“turnover mass” hint that the stellar mass of a galaxy 
plays a crucial role in quenching, possibly related to the 
“mass quenching” discussed in Peng et al. (2010). Fur¬ 
ther study is needed to determine the physical mecha¬ 
nism^) behind quenching. 

5.2. Increasing SFR with Redshift 

From our fits, we find strong evolution in So, which 
parameterizes the overall scaling of the SFR/M * main- 
sequence with redshift. 15 The scaling of the main se¬ 
quence has been found in the literature to evolve as 
(1 + z ) n , with the exponent n varying from 2.2 < n < 5 
(Erb et al. 2006; Daddi et al. 2007; Damen et al. 2009; 
Dunne et al. 2009; Pannella et al. 2009; Karim et al. 
2011). The value of n = 4.12±0.10 we measure is among 
the steeper slopes seen in the literature. 

It’s thought that the redshift evolution of the main- 
sequence normalization is due, at least in part, to in¬ 
creasing gas content in galaxies at earlier cosmic times. 
However, measuring the gas content in galaxies can be 
difficult, especially in high redshift systems. Molecu¬ 
lar hydrogen is notoriously difficult to detect, so many 
surveys instead probe the rotational transitions of CO 
and use locally calibrated CO-to-H2 conversion factors, 
although this conversion factor may differ in starburst 
galaxies (Tacconi et al. 2008; Magdis et al. 2011; Mag- 
nclli et al. 2012). Another method to estimate gas con¬ 
tent is to measure the dust mass from far-infrared or 
submillimeter photometry and convert to gas masses us¬ 
ing an assumed gas-to-dust ratio (e.g. Magdis et al. 2012; 
Santini et al. 2014; Scoville et al. 2014). Magdis et al. 
(2012) find gas fraction evolves as (1 + z) 2S , while re¬ 
cent ALMA observations suggest a steeper evolution of 
(1T z) 5 ' 9 (Scoville et al. 2014). Zahid et al. (2014) study 
the mass metallicity relationship and infer a much shal¬ 
lower evolution of gas mass M g oc (1 + z) 135 . Future 
studies will be needed to determine if the evolving nor¬ 
malization of the main-sequence can be explained simply 
by an increasing gas supply in galaxies, or if other expla¬ 
nations such as increased merger rates or increased star 

15 Alternatively, measuring the redshift evolution of the SFR at 
a constant characteristic mass provides similar results 
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formation efficiency are necessary to fully explain the ob¬ 
served evolution. 

6. CONCLUSIONS 

Using new far-infrared data from Herschel , we com¬ 
pare direct measurements of unobscured and obscured 
SFR with various SFR indicators that estimate the ob¬ 
scured SFR from data at shorter wavelengths (usually 
in optical or UV), and find that the NRK method of 
Arnouts et al. (2013) provides the most consistent esti¬ 
mate of the far-infrared derived SFR. By combining the 
SFRs from Herschel , Spitzer , and NRK, we analyze the 
relationship between SFR and M* (commonly referred 
to as the “star-forming main-sequence”) in 62,521 star¬ 
forming galaxies at 2 < 1.3 in the COSMOS field. From 
our new analysis we find: 

• The relationship between SFR and stellar mass 
does not follow a simple power-law, but flattens 
to near-constant SFRs at high stellar masses. The 
shape of the main sequence is roughly constant for 
all redshifts z < 1.3. 

• The scaling of the entire star-forming main se¬ 
quence rises with redshift as (1 + 2) 4124=01 °. 

• The characteristic turnover mass lies at Mo ~ 
10 10 Mq, with possible evolution toward higher 
turnover masses at high redshift. 

• The slope of the low-mass power-law lies between 
7 = 0.9-1.3, with possible weak evolution toward 
steeper slopes at higher redshift. 

• A broken power-law fit to galaxies below and above 
the turnover mass results in SFR oc jVf0.88±o.06 
below the turnover mass and SFR oc M°' 2 ' ±0 - 04 
above the turnover mass. 

Our analysis suggests that star-forming galaxies cannot 
be described by a single power-law relationship between 
SFR and M*, as had been suggested in many previous 
studies. Because of the strong effects of dust, direct ob¬ 
servations in the FIR are crucial for studying the entire 
population of star-forming galaxies. In future work we 
will explore possible causes of the turnover in the main 
sequence by studying detailed morphology and examin¬ 
ing possible feedback mechanisms, and we will extend 
our analysis to higher redshifts. 
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Table 1 


Redshift Range 

< z > 

log(Mnmit) 
log (Mo) 

So 

log(M e /yr) 

Mo 

log(M Q ) 

7 

Reduced \ 2 

0.25-0.46 

0.36 

8.50 

0.80 ± 0.019 

10.03 ± 0.042 

0.92 ± 0.017 

1.74 

0.46-0.63 

0.55 

9.00 

0.99 ± 0.015 

9.82 ± 0.031 

1.13 ± 0.033 

1.52 

0.63-0.78 

0.70 

9.00 

1.23 ± 0.016 

9.93 ± 0.031 

1.11 ± 0.025 

1.48 

0.78-0.93 

0.85 

9.30 

1.35 ± 0.014 

9.96 ± 0.025 

1.28 ± 0.034 

1.84 

0.93-1.11 

0.99 

9.30 

1.53 ± 0.017 

10.10 ± 0.029 

1.26 ± 0.032 

0.62 

1.11 1.30 

1.19 

9.30 

1.72 ± 0.024 

10.31 ± 0.043 

1.07 ± 0.028 

1.24 


Note. — Parameters of the best fit model to the star-forming main sequence. The full sample of 62,521 star-forming 
galaxies is split into six equally populated bins, with each bin containing ~ 17745 galaxies. Within each redshift bin, the 
galaxies are split into 30 equally populated bins of stellar mass. The median SFR in each mass bin is calculated and then 


fit to *5 = <Sq — log 


1 + 




( io- 

V 10^0 


, where S = log(SFR) and A4 = log(M *). Table columns are as follows: (1) Redshift 


range of bin; (2) Median Redshift; (3) Stellar Mass Limit of redshift bin; (4) So, the maximum value of S; (5) Turnover 
Mass; (6) Low-mass power-law slope; (7) Reduced \ 2 °f fit- 
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Figure 1 . Comparison of total SFR determined from combining direct measurements of FIR (Herschel , Lee et al. 2013) and UV (GALEX, 
Zamojski et al. 2007) with various SFR indicators using observations at shorter wavelengths. The different SFR indicators are: [top left] 
Spitzer 24 pm (Le Floc’h et al. 2009) + GALEX; [top right] multi-wavelength SED fits (Ilbert et al. 2013); [bottom left] NRK (0 < z < 1.3, 
Arnouts et al. 2013); and [bottom right] BzK (1.4 < 2 < 2.5, Daddi et al. 2007). In each panel, black contours give the density and 
concentration of sources, with extreme outliers plotted as gray circles. We bin the data in 20 equally populated bins (except BzK, which 
has 8 bins) and find the median SF Rindicator in each bin. Errors on the median points are measured using a bootstrapping technique and 
are plotted when larger than the size of the symbol. At the top of each panel is the Pearson correlation coefficient (with a value of +1 
indicating strong positive correlation and 0 indicating no correlation) and the typical difference between each particular SFR indicator and 
the Herschel -derived SFR. 
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Figure 2. Median difference in SFR estimated from infrared vs. optical/UV indicators, as a function of Total SFR. Sources that were 
detected by both SFR indicators were split in 15 equally populated SFR^otai bins, and the median value is plotted, with error bars 
representing bootstrapped errors. We see that all indicators begin to deviate significantly at log(SFR) > 1.5. (Top) Comparison of 
Spitzer 24 pm + GALEX (blue), multi-wavelength SED fits (green), NRK (orange, 0 < z < 1.3), and BzK (red, 1A < z < 2.5) to total 
SFR as derived from Herschel. (Bottom) Comparison of multi-wavelength SED fits (green), NRK (orange, 0 < z < 1.3), and BzK (red, 
1.4 < z < 2.5) to total SFR as derived from Spitzer 24 pm. 
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Figure 3. Relative fraction of sources with SFRs determined from NRK (dotted line), Spitzer 24 pm (dashed line), and Herschel (solid 
line). Different colored lines represent percentages in different redshift bins, with bluer colors representing low redshifts and redder colors 
representing high redshifts. 
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Figure 4. Contour density plot of star-forming galaxies in the COSMOS field. We remove all galaxies classified as “quiescent” (unless 
they were detected in the infrared) and combine all star-forming galaxies, regardless of the specific SFR indicator used. To display the 
density of sources in the SFR/M* plane, each redshift slice was made into a grid of 51 X 51 bins, and the number of sources in each bin was 
calculated. Black contours show the density of galaxies, with contour levels set at 1/2 of the standard deviation of the number of sources 
in each bin. Colored vertical bars represent the median SFR in mass bins of width 0.3 dex and display the overall trend of the SFR/M * 
relationship. Histograms of matching color display the distribution of SFR in each mass bin along the sides of each plot. Main-sequence 
relationships from Karim et al. (2011, green dots) and Whitaker et al. (2014, blue line) are plotted for comparison. 
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Figure 5. Median SFR in 6 equally populated redshift bins that have been split into 30 equally populated stellar mass bins. Errors on the 


median are calculated from bootstrapping. Solid lines represent the best-fit curve to the model S = So — log 
dashed line represents the stellar mass limit below which NRK has not been well-calibrated. 
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Figure 6. Top: Redshift evolution of the best-fit parameter values for (A) maximum of log(SFR), So', (B) turnover mass, Mo', and (C) 
power-law slope, 7 . Different color dots represent parameter values in different redshift bins with 1-cr error bars, and the dotted line is the 
best-fit linear fit to the data. Bottom: 95% confidence error ellipses displaying the covariance between (D) So and Mo', (F) So and 7 ; 
and (F) 7 and Mo, with different colors once again representing different redshift bins. We see moderate covariance in 7 and Mo, but 
little covariance between the other pairs. 



M(r)-M(J) 

Figure 7. NUV — r - *" vs. — J plot of a mass-selected sample of galaxies in COSMOS at 0.2 < z < 1.3. Black contours represent the 
full sample of galaxies, while orange circles highlight galaxies that are IR-detected, either with Herschel or Spitzer 24 ^m. The red line, 
from 113, divides the sample into “star-forming” and “quiescent” galaxies. The majority of IR-detected galaxies are properly classified as 
“star-forming”, but there is a significant population 7%) that are misclassified as “quiescent”. 
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